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ABSTRACT 

We have carried out VLA Zeeman observations of H I absorption lines toward the 
H II region in the Ml 7 giant molecular cloud complex. The resulting maps have 60"x45" 
spatial resolution and 0.64 km s _1 velocity separation. The Hi absorption lines toward 
M17 show between 5 and 8 distinct velocity components which vary spatially in a 
complex manner across the source. We explore possible physical connections between 
these components and the Ml 7 region based on calculations of H I column densities, line 
of sight magnetic field strengths, as well as comparisons with a wide array of previous 
optical, infrared, and radio observations. 

In particular, an H I component at the same velocity as the southwestern molecular 
cloud (M17 SW) ~ 20 km s _1 seems to originate from the edge-on interface between the 
H II region and M17 SW, in un-shocked PDR gas. We have detected a steep enhance- 
ment in the 20 km s _1 Hi column density and line of sight magnetic field strengths 
(Bi os ) toward this boundary. A lower limit for the peak 20 km s _1 H I column density is 
NhiIT s > 5.6 x 10 19 cm" 2 / K whil e the peak Bi os is ~ —450 /iG. In addition, blended 
components at velocities of 11—17 km s _1 appear to originate from shocked gas in the 
PDR between the Hll region and an extension of M17 SW, which partially obscures 
the southern bar of the Hll region. The peak Nhi/T s and Bi os for this component are 
> 4.4 x 10 19 cm~ 2 /K and ~ +550 fiG, respectively. Comparison of the peak magnetic 
fields detected toward M17 with virial equilibrium calculations suggest that ~ 1/3 of 
M17 SWs total support comes from its static magnetic field and the other 2/3 from its 
turbulent kinetic energy which includes support from Alfven waves. 



Subject headings: Hll regions — ISM:clouds — ISM:individual (M17) — ISM:magnetic 
fields — radio lines: ISM 
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1. INTRODUCTION 

In recent years it has become increasingly 
clear that magnetic fields play an important role 
in the process of star formation. (See, for ex- 
ample, Mouschovias &; Spitzer 1976; Heiles et 
al. 1993; Mckee et al. 1993). Unfortunately, 
the experimental techniques for measuring the 
strength and direction of magnetic fields are few 
and these methods are observationally challeng- 
ing. One such technique uses the Zeeman effect 
in 21 cm H I absorption lines toward galactic H n 
regions. Very Large Array (VLA) observations 
of the Zeeman effect in this line yield maps of 
the line-of-sight magnetic field strengths (Bi os ). 
These maps can then be compared to maps of the 
distribution and kinematics of ionized, atomic, 
and molecular gas to study the role of magnetic 
fields in star forming regions. Moreover, esti- 
mates of the masses, column densities, and vol- 
ume densities of interstellar material associated 
with the star-forming regions can be compared to 
field strengths measured via the Zeeman effect. 
These comparisons are needed to determine the 
energetic importance of magnetic fields in star- 
forming regions. For studies of this type, the M17 
H II region— molecular cloud complex is ideal be- 
cause it has been extensively studied at many 
wavelengths. 

The M17 Hn region is one of the strongest 
thermal radio sources in our galaxy. Its dis- 
tance has been estimated through photometric 
and kinematic means to be ~ 2.2 kpc (Chini, 
Elsasser, & Neckel 1980; Reifenstein et al. 1970; 
Wilson et al. 1970). There are at least 100 stars 
in the M17 Hil region with one 04 V (Klein- 
mann's star) and three 05 V stars providing the 
bulk of the ionizing radiation (Felli, Churchwell, 
& Massi 1984 and references therein; see Fig. 1). 
The H II region is made up of two distinct bar like 
structures which are ~ 5.7 pc long and ~ 1.1 to 
1.5 pc wide. The southern bar has nearly twice 
the peak radio continuum flux as the northern 
bar, and is located just to the east of the M17 SW 



molecular cloud core. 

The near absence of optical radiation from the 
southern bar indicates that this region suffers 
much more optical extinction than the northern 
bar (Gull & Balick 1974, Fig 1). Dickel (1968) 
and Gately et al. (1979) estimate A v ~ 1-2 mag 
toward the northern bar. Estimates of the vi- 
sual extinction toward the southern bar range 
from A v ~ 10 near the radio continuum peak 
to A v ~ 200 further west toward the core of the 
adjacent molecular cloud (M17 SW) (Beetz et 
al. 1976; Felli et al. 1984; Gatley et al. 1979; 
Thronson & Lada 1983). The low extinction to- 
ward the northern bar seems to indicate that the 
northern molecular cloud, first identified by Lada 
et al. (1976) (vlsr ~ 23 km s" 1 ), is behind the 
Hil region (Chrysostomou et al. 1992). 

The most interesting characteristic of M17 is 
that the interface between the southern bar of 
the Hn region and M17 SW is seen almost edge 
on (Fig. 1). From models of the density gradi- 
ent on the western edge of the H n region, Icke, 
Gatley, & Israel (1980) estimated that this in- 
terface lies at an angle of ~ 20° with respect to 
the line of sight. Although the transition from 
the Hil region to the molecular cloud is quite 
sharp, the interface region is too wide in pho- 
todissociation region (PDR) tracers such as [C n] 
158 /im unless it is clumpy (see eg., Stutzki et 
al. 1988). The extended [Cn] 158 /im emission 
observed in M17 SW by Stutzki et al. (1988) 
and Boreiko, Betz, & Zmuidzinas (1990) requires 
912-2000A UV photons to persist into the molec- 
ular cloud at least an order of magnitude fur- 
ther than the predicted UV absorption length 
scale. That is, atoms are being photoionized 
and molecules photodissociated further into the 
molecular cloud than one would expect for a ho- 
mogeneous medium. (See Tielens & Hollenbach 
1985 for homogeneous model calculations.) A 
possible explanation for this paradox is that the 
interface region is clumpy. 

Evidence that M17 SW is clumpy has been 
observed with high angular resolution at many 
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RIGHT ASCENSION (1950) 

Fig. 1. — Map of the 1, 5, 10, 20, and 30 Jy beam -1 21 cm continuum contours (present work, (thick 
lines) superimposed on an optical image from Gull & Balick and simplified 12 CO contours (thin lines) 
from Thronson & Lada (1983). The four stars which are primarily responsible for ionizing the Hll region 
are indicated with star symbols, with the largest representing Kleinmann's star (Felli, Churchwell, & 
Massi 1984) . 
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wavelengths. For example, Felli et al. (1984), 
found dumpiness in their ~ 10" resolution radio 
continuum observations of the M17 Hn region. 
They suggest this dumpiness originates from 
high density neutral clumps surrounded by dense 
ionized envelopes. These clumps may be the rem- 
nants of clumps from the parent molecular cloud 
which have been overtaken by the H n region ion- 
ization front. Massi, Churchwell, & Felli (1988) 
observed seven clumps in NH3 emission with 6" 
resolution toward the northwestern portion of the 
Hn region— M17 SW interface. The average den- 
sity and size of these clumps are n# 2 ~ 10 7 cm -3 
and 0.05 pc. Stutzki & Glisten (1990) modeled 
their ~ 13" resolution C 18 (2 — > 1) observa- 
tion (with a larger field of view than the NH3 
observation) with ~ 180 clumps with average 
densities of n# 2 ~ 10 5-6 cm -3 and diameters 
of ~ 0.1 pc. Hobson (1992) also saw evidence 
for dumpiness in HCO + and HCN molecules at 
~ 19" resolution, and Wang et al. (1993) found 
clumps consistent with the C 18 clumps in ~ 18" 
to ~ 34" resolution observations of CS and C 34 S 
molecules. 

The identification of clumps in M17 SW is sig- 
nificant for several reasons. For example, calcu- 
lations of density are more complex due to its 
dependence on the optical depth of the observed 
species, and the assumed clump filling factor. As 
discussed previously, the lengthscale of the PDR 
also depends on the dumpiness of the region. 
Therefore, the regions and morphology of pho- 
todissociated atomic gas are difficult to predict. 
Also the random motion of such clumps can in- 
fluence observed linewidths and the sizescale of 
turbulence within the molecular cloud. These 
considerations play an important role in our anal- 
ysis of the Hi gas and the line of sight magnetic 
fields calculated from it. 

In this paper we report the characteristics of 
the 21 cm continuum (§3.1), Hi optical depths 
and column densities (§3.2), and the line of sight 
magnetic field strengths (§3.3) derived from our 
VLA Hi Zeeman effect observation. We dis- 



cuss possible origins for the most prominent Hi 
absorption components in §4.1. In addition, 
we present comparisons of our 20 km s -1 B[ os 
map with previous linear polarimetry observa- 
tions in §4.2, and we also compare expected mag- 
netic field strengths from virial equilibrium argu- 
ments to our observed line of sight magnetic field 
strengths in §4.3. Our findings are summarized 
in §5. 

2. OBSERVATIONS 

Our Hi absorption data were obtained with 
the DnC configuration of the VLA []. Key pa- 
rameters from this observation can be found in 
Table 1. We observed both senses of circular po- 
larization simultaneously. Since the Zeeman ef- 
fect is very sensitive to small variations in the 
bandpass, we switched the sense of circular polar- 
ization passing through each telescope's IF sys- 
tem every 10 minutes with a front-end transfer 
switch. In addition, we observed each of the cal- 
ibration sources at frequencies shifted by ± 1.2 
MHz from the Hi rest frequency to avoid con- 
tamination from Galactic Hi emission at veloci- 
ties near those of M17. 

The calibration, map making, cleaning, and 
calculation of optical depths were all carried out 
using the AIPS (Astronomical Image Processing 
System) package of the NRAO. The right (RCP) 
and left (LCP) circular polarization data were 
separately calibrated and mapped using natu- 
ral weighting. These two data cubes were then 
combined to create Stokes I (RCP + LCP) and 
Stokes V (RCP - LCP) cubes. The I and V 
cubes were CLEANed using the AIPS task SD- 
CLN down to 70 mJy beam -1 . Bandpass correc- 
tion was applied only to the I data since band- 
pass effects subtract out to first order in the 
V data. Subsequent antenna leakage correction 
and magnetic field derivations were carried out 

lr The National Radio Astronomy Observatory is operated 
by Associated Universities, Inc., under contract with the 
National Science Foundation. 
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using the MIRIAD (Multichannel Image Recon- 
struction Image Analysis and Display) processing 
package from BIMA. 

3. RESULTS 

3.1. The 21 cm Continuum 

Figure 1 shows our 21 cm continuum contours 
with an optical image from Gull k, Balick (1974) 
and simplified 12 CO contours from Thronson & 
Lada (1983) superimposed. The total continuum 
flux at 21 cm is ~ 800 Jy calculated inside the 
1 Jy beam -1 (250 K) contour level (lowest con- 
tour in Fig. 1). This value is almost twice that 
obtained by Lockhart & Goss (1978) with a 2' 
beam at the Owens Valley interferometer. In ad- 
dition, we estimate the total flux of the South- 
ern bar is 384 Jy while that of the Northern bar 
is 240 Jy (calculated for points inside the 5 Jy 
beam -1 contour level). Our total flux inside the 
5 Jy beam -1 contour (624 Jy) is similar to the 
644 Jy reported by Lada et al. (1976) with the 
Haystack 120 ft telescope and 620 Jy observed by 
Lobert & Goss (1978) with the Fleurs synthesis 
telescope at 21 cm. 

3.2. Atomic Hydrogen Optical Depths and 
Column Densities 

M17 H I optical depths were calculated assum- 
ing that the spin temperature T s is much less 
than the background continuum temperature T c 
of the H II region. In addition, optical depth cal- 
culations were restricted to positions where the 
continuum power is greater than 1 Jy beam -1 
(250 K) and the signal-to-noise ratio is higher 
than three. The complex nature of the absorp- 
tion lines in M17 can be seen from the opti- 
cal depth spectra displayed in Figure 2. The 
strengths of some components vary spatially over 
the source but remain distinct, while other com- 
ponents seem limited to particular regions. This 
complexity makes Gaussian fitting very difficult. 
Our profiles agree qualitatively with the Gaus- 
sian fits reported by Lockhart & Goss (1978) 



made with a 2' beam and 0.85 km s _1 velocity 
resolution. They fit their M17 Hi optical depth 
profiles with eight Gaussians having center ve- 
locities ranging from 4.2 to 27.5 km s _1 . In the 
following analysis we concentrate on the Hi ve- 
locity component near 20 km s -1 and the blended 
components between 11 — 17 km s -1 due to their 
high optical depths and their coincidence in ve- 
locity with other species observed in M17 SW 
(see §4.1 for details). 

Figure 3 displays the morphology of the Hi 
optical depths toward the M17 Hn region as a 
function of velocity (every two adjacent channels 
were averaged). The steep increase of the Hi 
optical depth toward the Hn region— molecular 
cloud interface is clearly apparent. This effect 
is particularly noticeable near 20 km s _1 . Al- 
though the Hi optical depths become saturated 
at the western boundary of the Hn region (par- 
ticularly in the northwest), they reach values at 
least as high as 5. Optical depth maps in the 
11 — 17 km s -1 range have their highest values 
further to the east and in the northern part of 
the source. 

The column density of Hi toward M17 was 
found using the relationship 

N HI = 1.823 x 10 18 T S J r u du cm" 2 , (1) 

where r v is the optical depth per unit frequency. 
Values for Nhi/T s summed across the entire Hi 
velocity range, the 20 km s -1 component (from 
17 to 24.5 km s _1 ) and the 11—17 km s" 1 blended 
component are displayed in Figs. 4a, 4b, and 4c, 
respectively. At positions along the northwest 
and southwest portions of the source, these maps 
represent lower limits to Nhi/T s owing to satura- 
tion effects. Note the increase in the Hi column 
density toward the Hn region— molecular cloud 
interface on the southwestern side of the maps. 
There is also spatial agreement between the Hi 
column density concentration prominent in the 
total N HI /T S map (Fig. 4a) near 18 h 17 m 30.0 s , 
— 16°13'05" and the northern condensation seen 
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Fig. 2. — Hi Optical depth profiles across the source after convolution with a 2' beam (each profile is 
independent). The 1 Jy beam -1 21 cm continuum contour, and 2' beam are shown for reference . 
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Fig. 3. — Velocity channel maps of the H I optical depth morphology (after averaging every two channels) 
from 6.0 to 25.3 km s _1 (greyscale), with contours at r =1, 2, 3, 4, and 5. Notice the increase in Hi 
optical depths as you approach M17 SW in the 24.0 to 17.6 km s _1 velocity range . 
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Fig. 4. — Maps of Hi column density divided by spin temperature (Nhi/T s ) calculated using equation 
(1). For each map the outermost black contour is the 1 Jy beam -1 21 cm continuum level, and the 
greyscale and white contours represent Nhi/T s x 10 18 cm~ 2 /K. (a) Data summed across the entire Hi 
line. White contour levels are in increments of 10 x 10 18 cm _2 /K from 20 to 100 xlO 18 cm _2 /K, (b) Data 
for the 20 km s _1 component summed from 17 to 24.5 km s _1 . White contour levels are in increments of 
5 xlO 18 cm _2 /K from 5 to 55 xlO 18 cm~ 2 /K, and (c) Data for the 11 to 17 km s" 1 component summed 
from 10.5 to 16.9 km s _1 . White contour levels are in increments of 5 xlO 18 cm~ 2 /K from 5 to 40 xlO 18 
cm- 2 /K . 

7 



in molecular gas (see Wang et al. 1993 and 
Fig. 8). The total Nhi/T s at this position is 
8.3 xlO 19 cm- 2 /K. 

Aside from line saturation, the main source of 
uncertainty in calculating Hi column densities 
is in the value assumed for T s . Upper limits to 
T s come from two sources. For one, the H I line 
widths in some cases are quite narrow. Toward 
the southern end of the molecular ridge defined 
by CS emission (Wang et al. 1993) they are as 
little as 2.8 km s _1 for the 20 km s _1 component 
(Figs. 8, 9). This width implies a maximum tem- 
perature of 180 K at this position. Also, along 
the periphery of the continuum source, Hi ab- 
sorption lines exist where the continuum bright- 
ness is as little as 250 K (1 Jy beam -1 ). Evi- 
dently, T s is less than about 200 K in many direc- 
tions toward M17. At the position illustrated in 
Fig. 9, the CS line has a width of about 2 km s -1 . 
If H I and CS are mixed and at a common tem- 
perature, then T s «100 K. In fact, modeling of 
IRAS 12 fim infrared data byHobson & Ward- 
Thompson (1994) require a hot dust component 
of ~108 K to explain the overall 12 /xm emission 
from M17. This hot dust can be linked to the 
PDR by the fact that warm emission at 12 fj,m is 
likely to arise from a population of polycyclic aro- 
matic hydrocarbons (PAHs) (Desert, Boulanger 
& Puget 1990), and PAH emission from PDR re- 
gions has been predicted by many authors (e.g. 
Hollenbach 1993). 

Another indication of temperature in the atomic 
gas comes from the observations of [C i] 370 fim 
by Genzel et al. (1988). They estimate a tem- 
perature of 50 K for [Cl] 370 /im gas near 20 
km s _1 . Molecular temperature determinations 
may also be relevant if some of the H I is mixed 
with molecular gas (§4.2). Bergin et al. (1994) 
estimate Tk = 50 K from observations of low- J 
CO transitions, while Harris et al. (1987) show 
that high J CO lines (7 — > 6, 14 — > 13) arise from 
PDR gas at about 250 K. Meixner et al. (1992) 
found that the spatial extent and intensities of a 
wide array of far-infrared and submillimeter cool- 



ing lines could be adequately modeled by a three 
component gas consisting of very dense clumps 
n ~ 5 x 10 5 cm -3 , intermediate density inter- 
clump gas n~5x 10 3 cm -3 , and a tenuous ex- 
tended halo n ~ 3 x 10 2 cm~ 3 . They obtained the 
best agreement with observations by assuming an 
interclump gas kinetic temperature of ~ 200 K 
where it is also likely that dissociated H I gas can 
be found. Almost certainly, T s in the Hi gas is 
quite variable, with values in the range 50 - 200 
K most common. 

As described above, it is likely that no one 
spin temperature is applicable to all of the Hi 
gas components seen in M17. With a spin tem- 
perature of 50 K, a lower limit to the peak 
total column density toward M17 is Nhi = 
5.1 x 10 21 cm~ 2 . Spin temperatures between 
50 K and 200 K yield peak column densities 
of 2.8 - 11.0 x 10 21 cm~ 2 for the 20 km s" 1 
component and 2.2 — 9.0 x 10 21 cm -2 for the 
11 — 17 km s _1 blended components. We can use 
such estimates of the Hi column density to es- 
timate what fraction of the total hydrogen col- 
umn density arises from atomic gas. For exam- 
ple, Felli et al. (1984) used the optical depth of 
the silicon absorption feature at 10 /im to esti- 
mate that A u =30 toward the ultra compact Hn 
region M17-UC1 (18 h 17 m 32 s , -16°13'00"). Us- 
ing N Htot = 2 x 10 21 A V ( Bertoldi& Mckee 
1992), we estimate that the total column den- 
sity of hydrogen is Nh tot ~ 6 x 10 22 cm -2 . 
Our total Hi column density toward M17-UC1 
is N H i > 3.7 x 10 21 cm" 2 (T s =50 K). Unless 
saturation effects have resulted in a severe un- 
derestimate of Nhi, much of the gas along the 
line of sight to M17— UC1 is in the form of H2. 

3.3. Magnetic Field Strengths From the 
Zeeman Effect 

The Hi Zeeman effect is only sensitive to 
the line of sight magnetic field (Bi os ) in most 
astrophysical situations since the Zeeman split- 
ting is a tiny fraction of the line width. Since 
Voc dl/du, values for Bi os were obtained by fit- 
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ting the derivative of the I profile to the V profile 
at each pixel with a least squares fitting routine 
(for details see Roberts et al. 1993, Crutcher 
et al. 1996). Due to the complexity of the Hi 
spectra, the number of channels fit for each ve- 
locity component is small. To obtain more re- 
alistic error estimates, forty channels outside of 
the line region were included in each fit. The 
resulting Bi os maps were subsequently masked 
where the continuum power is less than 1 Jy 
beam -1 and when Bi os /a(Bi os ) < 3. The quan- 
tity Bi os /a(Bi os ) will subsequently be denoted 
S/N B . 

3.3.1. Magnetic Fields Derived for the 20 km s~ 1 
Component 

Observations of M17 have been made for many 
astrophysically relevant molecular and atomic 
transitions. Prom these data it is known that 
M17 SW has a LSR velocity of 19 to 21 km s -1 
which corresponds to the velocity of one of our 
strongest Hi absorption components. Figure 5 
shows an Hi optical depth profile and several 
molecular emission profiles from Stutzki et al. 
(1988) for a position at the H II region— molecular 
cloud interface. Note the similarity in the center 
velocities (~ 20 km s -1 ) of each species. This 
velocity coincidence strongly suggests that the 
20 km s -1 Hi absorption component is closely 
associated with the molecular cloud, most likely 
arising in the PDR at the Hll region-molecular 
cloud interface. Further evidence for this associ- 
ation comes from the steep rise in the 20 km s -1 
H I optical depth toward the interface region (Fig. 
4b). Therefore, values of B[ os derived from the 
H I Zeeman effect in this component are of direct 
relevance to the energetics of the molecular cloud 
and its immediate surroundings (see§4.1.2). 

A significant magnetic field exists in the 20 
km s _1 component over much of the central re- 
gion of the southern bar. The velocity range of 
the fit for this component is 17.0 to 24.6 km s -1 . 
A distinct Hi component at ~ 23 km s" 1 is 
also visible toward the northwestern corner of 



the source. A component at this velocity was 
also detected at this velocity by Lada & Chais- 
son (1975) in H2CO. However, detailed inspec- 
tion of the H I data cube shows that the channel 
range used to fit the 20 km s _1 component was 
such that no magnetic field was detected in re- 
gions where the 23 km s^ 1 component is distinct. 
Therefore, the Bi os detected in the velocity range 
17 to 24.6 km s _1 is only representative of the 
field in the 20 km s _1 Hi component. A map of 
Bi os for the 20 km s _1 component is shown in Fig. 
6, and sample fits at representative pixels are 
shown in Figs. 7a, 7b, and 7c. Negative values 
indicate that B[ os is toward the observer. Values 
of Bi os for this component (with sufficiently high 
S/Nb to pass our cutoff of 3) are centered near 
the continuum peak of the southern bar. The 
values for Bi os start out at ~ —100 fiG on the 
eastern side of the southern bar, pass through 
a shallow saddle point and then rapidly rise to 
~ —450 /uG on the western side of the map near 
the Hll region— molecular cloud interface. The 
maximum S/Nb obtained for this component is 
8.5. The apparently rapid change in direction of 
the Bi os along the narrow band in the northwest 
portion of the map may be real but is approach- 
ing the limit of our sensitivity and will require 
future observation to confirm. 

Note that the region of highest measured Bi os 
coincides closely with the ridge of CS emission 
observed by Wang et al. (1993; Fig. 8). Also, 
the 20 km s -1 Hi line in this region is very nar- 
row, and its width and center velocity correspond 
closely to those of the CS lines (Fig. 9). This ve- 
locity correspondence, like those of Fig. 5, pro- 
vide further evidence of a close association be- 
tween the 20 km s _1 Hi gas and the molecular 
gas. In addition, the narrowness of the lines in 
this region seem to exclude the possibility of the 
ridge being post-shock gas. The nature of the H 1 
gas at 20 km s -1 and its effect on the morphology 
of Bios will be discussed in §4.1.2. 
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Fig. 5. — An Hi optical depth profile (present work) and 13 CO, C 18 0, and NH3 emission profiles from 
Stutzki et al. 1988 toward the position 18 h 17 m 30.6 s , -16 o 14'00" (also known as position (-60,-30) in 
Stutzki et al.) . 
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Fig. 6. — Greyscale map of Bi os for the 20 km s _1 component with white contours at —400, —300, —200, 
— 100, and —50 [iG. The 21 cm 1, 5, 10, and 20 Jy beam -1 continuum contours are shown in black. The 
three white + symbols show the positions of the profiles shown in Figs. 7 a, b, c . 
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Fig. 7. — Representative fits at three positions for the 20 km s _1 component (17.0 to 24.6 km s -1 ). The 
upper panels show Stokes I profiles (dashed histogram), and the bottom panels show Stokes V profiles 
(solid histogram) with the fitted derivative of Stokes I shown as a smooth dashed curve. In the upper 
and lower panels, the velocity range used in the fit is denoted by solid or heavier lines. The value of B\ os 
fit for each position and its calculated error are given at the top of each plot, (a) Position 18 h 17 m 35 s , 
-16°14'00". (b) Position 18 h 17 m 30.6 s , -16°14'00", or position (-60, -30) from Stutzki et al. (1988). (c) 
Position 18 h 17 m 28.5 s , -16°14'16", or position (-90, -45) from Stutzki et al. (1988) . 
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RIGHT ASCENSION (.1950) 

Fig. 8. — Magnified greyscale of the Bi os for the 20 km s _1 component with black CS J = 2 — > 1 contours 
from Wang et al. (1993) superimposed. The white Bi os contours are at —400, —300, —200, —100, and 
—50 //G. Notice the spatial coincidence between the highest values of Bi os and the molecular ridge traced 
by CS emission . 
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Fig. 9. — An Hi optical depth profile (present work) and emission spectra from several CS transitions 
toward the molecular ridge (18 h 17 m 32.4 s , -16°15'06") from Wang et al. (1993). See also Fig.8 . 
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3. 3. 2. Magnetic Fields Derived for the Blended 
11— 17 km s _1 Components 

Strong, primarily positive magnetic fields were 
also detected in blended velocity components be- 
tween 11 to 17 km s _1 . A map of this magnetic 
field is shown in Fig. 10 and the fit to a few rep- 
resentative pixels can be seen in Figs. 11a, lib, 
and 11c. There are four main regions of signifi- 
cant Bi os for this blended component which ap- 
pear primarily in the northeastern, north central 
and western parts of the map. The center veloc- 
ity of this component (defined by the velocity at 
which the derivative changes sign) is somewhat 
different in each of these regions. They range 
from 14.4 km s _1 in the northeast (Fig. 11a), 
13.7 km s _1 in the central northern feature (Fig. 
lib), 14.0 km s _1 in the northwest (Fig. 11c) 
and 12.4 km s" 1 in the southwest. One interest- 
ing feature of the Bi os for this component is that 
the direction of the field changes from positive 
(away from the observer) to negative and back to 
positive across the map from east to west (this 
trend can be seen in the fits shown in Figs. 11 
a, b, c). Bi os in the northwestern feature reaches 
values as high as +550 fj,G with S/Nb as high as 
14. 

Evidence for molecular gas near this veloc- 
ity has been observed in a number of species. 
For example Greaves, White, & Williams (1992) 
reported a velocity component in this range in 
their 12 CO (3 -► 2) and 13 CO (3 -► 2) spec- 
tra. Their component shifts from 12 km s _1 in 
the southeast to 16 km s _1 further to the south- 
west of their field of view (the eastern portion of 
their maps overlaps the the western side of our 
maps). The positions where their low velocity 
component is centered at 12 km s _1 corresponds 
to the southwestern region of our magnetic field 
map where the Hi line center is 12.4 km s _1 . 
Rainey et al. (1987) and Stutzki et al. (1988) 
also identified emission in this velocity range in 
12 CO (3 -» 2) and 12 CO (2 -» 1), respectively. 
Massi et al. (1988) observed that two of their 
seven NH3 clumps near M17-UC1 have velocities 



of 15.5 km s 1 and Stutzki et al. (1988) observed 
31 C 18 clumps (out of 180) in the 11-17 km s" 1 
velocity range . 

In addition to molecular gas, there is also [C n] 
158 /j,m emission at a velocity of ~ 17.8 km s _1 
(Boreiko et al. 1990). [Cn] 158 fim is an im- 
portant tracer of PDR's and should exist in the 
same zone of a PDR as Hi gas. The precise spa- 
tial correlation between these two atomic species 
for Ml 7 is unknown. However, the low reso- 
lution (3/4) map of [Cll] 158 ^m obtained by 
Matsuhara et al. (1989) does seem to follow the 
same morphology as our summed Nffj/T s (see 
Fig. 4a). A profile from one pixel at the [Cn] 
158 /im peak from Boreiko et al. (1990; spa- 
tial resolution 43") and Hi are shown in Fig. 
12. The [Cn] 158 //m profile appears to have 
a component at a lower velocity than their low- 
est reported velocity of 17.8 km s _1 . The authors 
noted that there may have been some contamina- 
tion from emission in the reference beam at this 
position. The possible significance of a lower ve- 
locity [C 11] 158 /im component will be discussed 
again in §4.1.3. 

3.3.3. Magnetic Field Detections at ~ 3.4, 6.7, 
and 27 km s^ 1 

Line of sight magnetic field detections were 
also made for three other velocity components. 
A significant B[ os was measured in the velocity 
range —2.4 to 4.0 km s _1 with the line center 
at ~ 3.4 km s" 1 . This detection is interesting 
because the S/Nb in the southeastern region is 
as high as 10 and there is a rapid change in 
Bi os from negative to positive values along the 
northwestern border with M17 SW. This abrupt 
change in field direction corresponds to the be- 
havior of the Bi os measured in the 20 km s _1 
component in this same region. The line of sight 
fields in this component range from ~ —200 to 
+100 fjtG. 

Detections were also made in the velocity 
range from 4.1 to 10.5 km s _1 with a center ve- 
locity of ~ 6.7 km s _1 and the maximum 
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Fig. 10. — Greyscale map of B\ os for the 11—17 km s blended components with white contours at 
-250, -150, 100, 200, 300, 400 and 500 /j,G. The 21 cm continuum contours are shown in black at 1, 5, 
10, and 20 Jy beam -1 . The three + symbols show the positions of the profiles shown in Fig. 11 a,b,c . 
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Fig. 11. — Same as Fig. 7 a,b,c except the range of velocities used in the fit for B\ os was 11—17 
km s _1 . (a) Northeastern position 18 h 17 m 40.6 s , -16 o 10'15" (white + symbol in Fig. 10). (b) North 
central position 18 h 17 m 35.6 s , -16°10'51" (white + symbol in Fig. 10), and (c) Northwestern position 
18 h 17 m 28.7 s , -16°09'12" (black + symbol in Fig.10). Notice how the sign of the field changes as you 
move from east to west across the map . 
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S/Nb = 6.4. Reliable measurements of Bi os 
in this Hi component range from ~ —125 to 
+W0jj,G. The regions of significant Bi os mea- 
sured for this component do not particularly cor- 
respond to the morphology of the 21 cm contin- 
uum emission of M17. In fact the highest fields 
and concentration of H I in this velocity range are 
to the southeast. One interesting possibility is 
that there is a connection between Hi measured 
in this velocity range and the 7 km s _1 and the 
cold cloud observed by Riegel & Crutcher (1972) 
at 21 cm over a large region toward the galactic 
center including the direction of M17. Crutcher 
& Lien (1984) estimate that this cold cloud (~20 
K) must be located within 150 pc of the sun. 

A sparse Bi os was also detected in the range 
24.7 to 33.0 km s _1 with a center velocity of ~ 
27 km s _1 and maximum S/Nb=Q. The Bi os 
in this Hi component ranges from ~ —200 to 
+100 fiG. A component in this velocity range 
was also observed by Greaves et al. 1992 in 12 CO 
(3 -» 2) and 13 CO (3 -» 2) emission. The ratio 
of these two species' intensity seems to indicate 
the presence of a distinct cloud at this velocity. 

4. DISCUSSION 

4.1. Proposed Origin of Magnetic Field 
Components 

4-1.1. General Nature of the Mil Region 

Meixner et al. (1992) suggest that the coinci- 
dence of part of their [Oi] 63 fim map with the 
radio continuum emission indicates that hot ion- 
ized gas from the Hll region is carving a bowl 
into the molecular cloud to the southwest (Ml 7 
SW). The symmetry axis of the bowl lies in the 
plane of the sky, and the overlapping [O i] 63 ^m 
emission arises from a PDR on the front and/or 
back sides of the bowl. The bottom of the bowl 
(at the southwestern edge of the southern bar) 
is the portion of the Hll region-molecular cloud 
interface viewed edge— on. Further evidence of a 
bowl— like morphology is the high visual extinc- 
tion toward the southern bar which presumably 



arises from the front side of the bowl. Gull & Bal- 
ick (1974) also suggest a bowl morphology based 
on the double peaks separated by ~ 20 km s _1 in 
their H76a radio recombination line data. (Also, 
see Goudis Sz Meaburn 1976.) These two recom- 
bination line velocity components are thought to 
arise from streaming motions of newly-ionized 
gas at the front and back sides of the bowl. 

In many ways, this M17 interface model re- 
sembles the interface region between the Orion 
nebula and the molecular cloud behind it. For 
Orion, the interface lies in the plane of the sky 
with the ionizing stars directly in front of it. For 
Ml 7, the geometry is turned by 90°. The inter- 
face layer is perpendicular to the plane of the sky, 
with the ionizing stars offset to the northeast. 
The densest accumulation of ionized gas in M17 
lies along the "working surface" near the ioniza- 
tion front. Therefore, the brightest region in the 
M17 continuum map, the southern bar, does not 
coincide in the sky with the ionizing stars as it 
does for the Orion nebula. 

The dumpiness of the M17 interface has been 
observed in many molecular species (§1). It 
seems to closely resemble the scenario discussed 
by Bertoldi& Draine (1996), whereby dense clumps 
are slowly photodissociated by a nearby Hn re- 
gion. In this case the H II region is partially em- 
bedded in (and partially obscured by) a bowl 
of clumpy neutral gas. The Hi gas associated 
with Ml 7 can arise, in principle, anywhere in 
front of the source. Our data suggest that the 
20 km s _1 component originates primarily in the 
edge— on portion of the bottom of the molecu- 
lar bowl. In this model the 20 km s _1 Hi arises 
from molecular gas that been partially photodis- 
sociated but otherwise dynamically undisturbed 
(unshocked). On the other hand, the Hi data for 
the 11—17 km s^ 1 component supports the idea 
that it arises from photodissociated, shocked gas 
from the front side of the molecular bowl which 
is streaming toward us. The association of dif- 
ferent Hi absorption components with shocked 
and unshocked gas has also been suggested for 



16 



S106 (Roberts et al. 1995) and W3 (Roberts et 
al. 1993). Evidence supporting this picture for 
M17 is given in §4.1.2 and 4.1.3. 

4.1.2. The Nature of Hi Gas at 20 km s~ l 

The velocity range of the 20 km s _1 Hi com- 
ponent (17 to 24.6 km s _1 ) includes most of the 
molecular gas in the M17 region (e.g Rainey et al. 
1987 and §3.3.1). Therefore, Hi components in 
this range are likely to be closely associated with 
the bulk of the M17 molecular gas. Such an as- 
sociation is particularly clear toward the molec- 
ular ridge (see Fig. 8) where the center velocity 
of the 20 km s _1 Hi component closely matches 
that of a variety of molecular lines, including high 
density tracers such as CS (Figs. 5, 9, and 12). 
This Hi component is also relatively narrow (3- 
5 km s _1 ), and its center velocity changes little 
with position along the interface region. There- 
fore, it most likely originates in quiescent gas 
that is dynamically unaffected by the H II region. 
Much of this quiescent gas must lie in the edge- 
on Hn region— M17 SW interface region at the 
bottom of the bowl described by Meixner et al. 
(1992). Hence the concentration of 20 km s~ x Hi 
gas and high Bi os (~ —450 fiG) near the interface 
region (Figs. 3, 4b, 6, and 8). 

Further indication of the association between 
Hi and molecular gas comes from similarities 
in their spatial distribution near the molecular 
ridge. Fig. 13 (upper panel) shows the rela- 
tive intensity of several species along the NE- 
SW strip scan defined by Stutzki et al. (1988). 
The strip lies perpendicular to the interface re- 
gion and cuts across it (ionizing radiation propa- 
gates from left to right). As shown in the figure, 
the increase in the 20 km s" 1 N HI /T S and B los 
(lower panel) matches well the increase in 12 CO 
line strengths to the right of the Hll region. In 
fact, N HI /T S , 12 CO(2 -» 1), and 12 CO(7 -► 6) 
peak outside of the H n region as expected for an 
interface region viewed edge on. Note that the 
actual peak in Nhi/T s for the 20 km s" 1 com- 
ponent is not defined by our data since we can- 



not detect H I absorption beyond the continuum 
source. 

Notice that the integrated [C n] 158 /im emis- 
sion, also shown in Fig. 13, peaks inside (to the 
left of) the rise in the 20 km s" 1 N HI /T S . In 
a simple planar PDR model, [C n] 158 /im emis- 
sion coincides with H I and with molecular gas up 
to A v = 2 — 4 mag. (See review byHollenbach, 
1990.) However, the higher velocity resolution 
[Cn] 158 /im profiles of Boreiko et al. (1990) re- 
veal the presence of several [C n] 158 /im veloc- 
ity components. An examination of these pro- 
files suggests that the ~ 20 km s" 1 [C n] 158 /im 
component alone peaks further into the molecu- 
lar cloud than the integrated [C Ii] 158 /im emis- 
sion plotted in Fig. 13. Also, note that the rise 
in the 11-17 km s" 1 N HI /T S , plotted in Fig. 13 
(lower panel), peaks to the left of the 20 km s _1 
Nhi/T s in a similar manner to the integrated 
[Cll] 158 /jm emission. In fact, the contribution 
of lower velocity [C n] 158 /im emission to the in- 
tegrated emission plotted in Fig. 13 may explain 
why it seems to peak to the left of the 20 km s _1 
Nhi/T s . Apparently, Hi and [Cll] 158 /im gas is 
more broadly distributed along the strip in Fig. 
13 than is the gas at 20 km s _1 alone. This shows 
that absorbing H I gas lies on the near side of the 
bowl, not just at the bottom. 

The close association between H I gas near 20 
km s _1 and the molecular ridge region could arise 
in at least two ways. For one, the Hi may be 
well mixed with and a minor constituent of the 
molecular gas. In such a case, the Zeeman ef- 
fect near the ridge (Fig. 8) is directly represen- 
tative of magnetic fields in the molecular gas. 
A second possibility is that the absorbing Hi 
near 20 km s _1 is confined to a thin photodis- 
sociated shell about the molecular ridge or to 
thin shells surrounding a multitude of molecular 
clumps within the ridge. In this case, the H I Zee- 
man effect is also representative of the magnetic 
field in the molecular gas if the Hi is photodis- 
sociated but dynamically undisturbed H2. 

In either case, the H I Zeeman effect near the 
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Fig. 12. — An Hi optical depth profile (present work) and a [Cn] 158 spectra from Boreiko et al. 
(1990) for the position of the Stutzki et al. (1988) [Cn] 158 fjan "peak" 18 h 17 m 32.5 s , -16°13'42", or (-30, 
-15) . 
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Fig. 13. — Upper panel : Distribution of the relative intensity of 3.5 cm continuum, CO (7 — ► 6), CO 
(2 — > 1) and [C n] 158 /jm across the H II region— M17 SW interface following the NE-SW strip scan 
observed and plotted by Stutzki et al. (1988). Lower panel : Relative intensity of the 20 km s _1 Nhj/T s , 
20 km s _1 Bi os , and 11—17 km s _1 Nhi/T s for the same positions as the upper panel. The peak positions 
for the 20 km s _1 Nhi/T s and 20 km s _1 Bi os do not necessarily represent true maxima, due to the loss 
of continuum power as the strip is transversed to the west. Notice that the 11—17 km s _1 Nhi/T s 
distribution along this strip is similar to the distribution of [C II] 158 /im emission (east of the continuum 
power drop-off), while the 20 km s _1 Nhi/T s and 20 km s" 1 Bi os peak further west . 
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molecular ridge must arise in relatively dense gas 
if approximate equilibrium exists between the 
magnetic energy density and that associated with 
non-thermal motions in the gas. Such an equi- 
librium was suggested by Myers & Goodman 
(1988a, 1988b) on the basis of limited observa- 
tional data. In such Si CclSG, the gas density is 
given by 



B, 



avg 



cm 



where n e 



(2) 

is the aver- 



_0AAvnt. 
o eq is the proton density, B avg 
age field strength (/xG), and Avnt is the FWHM 
contribution to the line width from non-thermal 
motions (km s" 1 ; see also §4.3). As an exam- 
ple we apply this equation to the position of Fig. 
5 (~ 15" east of molecular ridge). At this po- 



4 km s 1 , and B, 



avg 



200 nG. 



sition, Av 

If T s = Tk = 50 K (as estimated by Bergin et 
al. 1994 from 12 CO observations), Av NT = 3.7 
km s _1 and n eq ^2x 10 4 cm~ 3 . This estimate 
of the local density in which the Hi Zeeman ef- 
fect arises is conservative since the actual field 
strength B avg almost certainly exceeds Bi os and 
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This value for n ea based on Bi 



is also in excellent agreement with the value of 
n eq «6x 10 4 cm -3 calculated toward M17 SW 
by Goldsmith, Bergin, & Lis (1997) using C 18 
column densities and average line widths. How- 
ever, at the position cited here, Nhi/T s in the 
20 km s _1 component is only 3 xlO 19 cm _2 K, so 
N HI = 1.5 x 10 21 cm" 2 for T s =50 K. A typi- 
cal scale length for H I absorption features in the 
plane of the sky is of order 2 pc (3' at 2.2 kpc). 
If the same scale length applies along the line of 
sight, then nui ~ 250 cm~ 3 . That is, if the Hi 
Zeeman effect arises in gas at a density of order 
then the Hi could be about only a 1% con- 
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stituent of largely molecular gas. This possibility 
exists because the Nh 2 estimated by Stutzki et 
al. (1988) and Goldsmith et al. (1997) near this 



direction 



10 



2.3 



cm 



is about 100 



times higher than Njji- Alternatively, Hi might 
not be mixed with H2 but confined to thin sheets 
or very dense, thin atomic envelopes surrounding 



the multitude of molecular clumps which have 
been observed in the interface region (see §1). 
Note that this would still require the Hi gas to 
occupy about 1% of the line of sight in the ab- 
sorbing region. These two cases are the same 
ones cited above to explain the close association 
of H I gas with molecular gas in the ridge. 

4.1.3. The Nature of Hi Gas at 11-17 km s _1 

The morphology and velocity structure of gas 
between 11—17 km s^ 1 suggest that it also lies 
very close to the Hn region on the near side of 
the bowl and that it has been shocked and ac- 
celerated away from the ionized gas. The 11—17 
km s _1 gas is primarily distributed in an arc that 
partially encircles the northern and southern bars 
of the Hn region (Figs. 3, 4c). This arc of Hi 
column density corresponds well with a similar 
arc structure apparent in 12 CO (3 — >2) emission 
mapped by Rainey et al. (1987) in the veloc- 
ity range 14—20 km s -1 . This morphology sug- 
gests a very close relationship between gas in the 
11—17 km s -1 velocity range and the Hn region. 
The gas is blue-shifted by 2-8 km s _1 relative 
to the Hn region center velocity of 19 km s _1 
(Joncas k, Roy 1986). Also, Hi in this veloc- 
ity range shows complex velocity structure with 
numerous components appearing at various posi- 
tions across the source (Figure 2, 11 a, b, c). This 
complicated velocity field and magnetic field re- 
versal (see §3.3.2) suggest that gas in the 11—17 
km s _1 range lies in an interface zone that has 
been dynamically affected by interactions with 
the Hn region, unlike the 20 km s _1 component 
discussed in §4.1.2. 

Overall, the association between atomic and 
molecular gas in this velocity range is much less 
clear than it is at 20 km s _1 near the west- 
ern interface region. This could in part be due 
to fact that few high resolution molecular ob- 
servations have mapped the northern regions of 
M17. However, judging from 12 CO J = 1 — >0, 
12 COJ = 3 ->2, and 13 CO observations (Lada 
1976; Rainey et al. 1987) which all show compo- 
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nents in this velocity range, the density of molec- 
ular material in this region is low compared to the 
western interface. Toward the northern bar in 
particular, where the 11 — 17 km s" 1 H I magnetic 
fields are high (see §3.3.2), there is not always a 
molecular counterpart in the 11 — 17 km s _1 ve- 
locity range. In fact, for the positions shown in 
Figs. 11a and lib there is no evidence of molecu- 
lar gas. Such a variable ratio of atomic to molec- 
ular gas might be expected if gas in this velocity 
range has been subjected to a widely variable de- 
gree of dissociation in the interface region. 

Although molecular column densities are low 
in the 11 — 17 km s" 1 range, Hi field strengths 
are very high, reaching values of over +500 //G 
(Fig. 10). Therefore, the assumption of equilib- 
rium between magnetic and non-thermal energies 
leads to conclusions similar to those for gas near 
20 km s" 1 (§4.1.2). That is, n eq > 10 4 cm" 3 , and 
riHi along the line of sight is about two orders 
of magnitude less than n eq . Yet the scarcity of 
molecular gas at positions of high field strength 
means that H I at these positions cannot exist as 
a minor constituent in a largely molecular region. 
If the H I alone is to exist at a density near 10 4 
cm -3 , it must reside in one or more thin layers 
that occupy only about 1% of the path length 
through the Hi absorbing region. That is, the 
absorbing H I gas is distributed in a very clumpy 
fashion along the line of sight. This conclusion is 
consistent with numerous other indications of a 
highly clumped medium in the M17 region (§1). 
Most likely, the absorbing H I gas is shocked, pho- 
todissociated gas lying in thin interface regions 
about ablating neutral clumps. The clumps, 
in turn, lie within mostly ionized gas, like the 
clumps modeled by Bertoldi & Draine (1996) or 
the clumps identified by Felli et al. (1984) in high 
resolution M17 radio continuum maps. In these 
clumps, the pressure of the ionized gas drives a 
shock front into the neutral medium. Therefore, 
a rough equivalence between the thermal energy 
density of the ionized gas and the magnetic en- 
ergy density in the adjacent compressed atomic 



region is expected. Evidently, this is the case. If 
the M17 ionized gas has a temperature of 8000 
K (Subrahmanyan & Goss 1996) and a density of 
10 4 cm" 3 (Felli et al. 1984), then its thermal en- 
ergy density is equivalent to the magnetic energy 
density in the thin Hi layers for B ~ 500 /uG. 

Several reasons exist to expect Hn region 
driven motions of a few km s _1 near M17. Ob- 
servationally, Rainey et al. (1987) fit their 13 CO 
data on molecular clumps near the southern bar 
to a geometric expansion model about Kleine- 
mann's star (Fig. 1). They obtain an expansion 
velocity of « 11 km s" 1 about this star. Ra- 
dio and optical recombination line observations 
often reveal split profiles shifted by about ±8 
km s _1 from the Hn region rest velocity (Gull 
& Balick 1974; Joncas & Roy 1986). Theo- 
retically, Bertoldi& Drain (1996) estimated the 
shock-induced acceleration expected in the Orion 
nebula PDR, obtaining a value of 3 km s _1 . This 
value is somewhat sensitive to the ratio of post 
to pre-shock density, taken by Bertoldi& Draine 
(1996) as 2. If the ratio for M17 is 4, as suggested 
by the ratio of the square of the 12 CO (7 — > 6) 
(5 km s" 1 ) line width (Harris et al. 1987) to the 
square of the CS (2 — > 1) line width (2.5 km s _1 , 
Wang et al. 1993), then the shock induced accel- 
eration is increased to 6 km s _1 . Also, Bertoldi) 
and Bertoldi h McKee (1990) discuss a rocket ef- 
fect upon ablating clumps of neutral gas exposed 
to ionizing radiation. They estimate that rocket 
acceleration should amount to a clump velocity of 
5—10 km s" 1 under typical circumstances. These 
estimates of shock speed, agree quite well with 
the 2-8 km s" 1 blueshift of the 11-17 km s" 1 
componets. 

4.2. Comparison of 20 km s 1 Bi os Field 
Morphology with Polarimetry Stud- 
ies 

A number of studies have been made of linear 
polarization toward M17. These include polar- 
ization at optical and near IR wavelengths from 
grain absorption and polarization in the far IR 
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from grain emission. (See references cited by 
Dotson 1996.) Such studies reveal the position 
angles of the magnetic field in the plane of the 
sky (B±), but they yield little if any information 
about field strengths. Therefore, they compli- 
ment Zeeman effect studies, assuming the linear 
and circular polarization arise in the same re- 
gions. In general, field directions inferred from 
grain absorption do not agree with those inferred 
from grain emission. Very likely, this disprepancy 
results from the different regions sampled by each 
technique, with grain emission data more likely 
representative of the field directions inside M17 
SW (Goodman et al. 1995). 

The study of polarized 100 /im dust emission 
by Dotson (1996) provides the most comprehen- 
sive available data set for comparison with Hi 
Zeeman effect results. In figure 14 we present a 
map of the position angle of the transverse mag- 
netic field -Bj_ from Dotson (1996) overlaid on 
our ~20 km s _1 Bi os grayscale. The resolution 
of the Dotson study (35") is comparable to that 
of the Hi data. Also, the 20 km s" 1 Hi compo- 
nent likely samples at least part of the molecular 
ridge (§4.1.2) as does the 100 micron emission. A 
principal conclusion of the Dotson study is that 
the M17 magnetic field shows a significant de- 
gree of spatial coherence across M17. The Hi 
Zeeman data for the 20 km s _1 component cor- 
roborates this result (Figs. 6, 8). Moreover, the 
eight-fold increase in Bi os from east to west to- 
ward the molecular ridge (Fig. 8) is matched by 
a decrease in fractional linear polarization from 
about 4% to 1%. Along this same line the field 
has a nearly constant position angle (predomi- 
nantly east- west). Since B[ os is perpendicular to 
both to the transverse magnetic field and the 100 
/im linear polarization vectors, one might expect 
to see a correlation between regions of maximum 
Bi os an d minimum 100 /im polarization. 

One obvious interpretation of the IR and H I 
Zeeman data is that the magnetic field lines curve 
into the line of sight as one approaches the molec- 
ular ridge from the east. In such a case, Bi os 



may be a close approximation to the total field 
strength in the region of the molecular ridge sam- 
pled by H I. Since the H I must lie in front of the 
Hn region, the field lines in this picture wrap 
around the Hn region on the front side of the 
bowl, becoming more aligned with the line of 
sight at the bottom of the bowl as they pass 
through the molecular ridge. This geometry is 
consistent with the inferences made by Dotson 
from the linear polarization data alone. The ori- 
gin of the field curvature may lie in the process of 
gravitational contraction that formed the molec- 
ular ridge, gathering an initially uniform field 
along the line of sight into an hourglass shape 
with the ridge at its waist. An hourglass shaped 
field structure was also suggested for the W3 re- 
gion by Roberts et al. (1993), although for this 
source, the axis of the hourglass is nearly in the 
plane of the sky rather than along the line of 
sight. Alternately, the present geometry of the 
M17 field may represent the bending effects of 
the Hn region as it expands into the molecular 
cloud deforming a field originally along the line 
of sight into a bowl-like shape. 

It should be noted that the occurrence of de- 
creasing linear polarization toward the the 100 
fj,m flux and column density maxima could re- 
sult from effects other than a true minimum in 
the transverse magnetic field component. For 
example, optical depth effects, dust properties, 
grain alignment efficiency, or some combination 
of these can cause a decrease in the linear polar- 
ization independent of the magnetic field. On the 
other hand, the magnetic field itself could simply 
be more tangled in higher density regions (Jones, 
Klebe, & Dickey 1992; Myers & Goodman 1991). 
In addition, the 760 /im polarimetry by Vallee & 
Bastien (1996) does not show a minimum in the 
linear polarization in this region of the molecular 
ridge. However, the 760 /tm polarization arises 
from cooler dust grains which could be spatially 
distinct along the line of sight from those emit- 
ting at 100 /im. Future polarimetry of this region 
at an intermediate wavelength from ISO and 
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Fig. 14. — Greyscale of the 20 km s Bi os overlaid with the 1, 5, 15 Jy beam contours from our 21 
cm continuum, and the 100 fim linear polarization vectors from Dotson (1996). The position angle of the 
100 fj,m polarization vectors indicates the direction of B±, which is perpendicular to the observed linear 
polarization direction . 
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higher resolution H I Zeeman results may help to 
clarify this issue. 

4.3. Magnetic Energy vs. Turbulent and 
Gravitational Energy 

The magnetic field in a cloud can be sepa- 
rately described in terms of a static component 
Bg and a time-dependent or wave component 
Bw- The static component connects the cloud 
to the external medium and also determines the 
total magnetic flux through the cloud, while the 
wave component is associated with MHD Alfven 
waves in the cloud (e.g. Arons & Max 1975). 
Although the static component of the field can 
only provide support to the cloud perpendic- 
ular to the field lines, the turbulent or wave 
component can provide 3-D support (Mckee & 
Zweibel 1995, also see references therein). If a 
cloud is threaded by = itR 2 Bs, it can be 
completely supported perpendicular to B$ if its 
mass does not exceed the magnetic critical mass 
M^^ cr it « 0.13G~ 1/2 $ B (Mouschovias & Spitzer 
1976). This criterion leads directly to the obser- 
vationally relevant relation 

B S ,crit ~ 5 x 10~ 21 N P nG, (3) 

where N p is the average proton column density of 
the cloud. If B s > B s , cr it (i-e. M < M$), then 
the cloud is supported by Bs, it is magnetically 
subcritical, and further evolution of the cloud 
perpendicular to the field can only occur via am- 
bipolar diffusion. If Bs < Bs )C rit, (i-e. M > M<j>), 
then Bg cannot fully support the cloud, the cloud 
is magnetically supercritical, and internal mo- 
tions must supply additional support if the cloud 
is stable. In this case, further evolution of such a 
cloud will be controlled in part by processes that 
create and dissipate internal motions. 

As discussed above, the second type of support 
provided by magnetic fields arises from Alfven 
waves generated in a turbulent velocity field. The 
fluctuating or wave component of the magnetic 
field Bw organizes supersonic (but sub-Alfvenic) 
motions in the gas so that highly dissipative 



shocks do not occur. Theoretical studies suggest 
that equipartition between the wave's kinetic 
and magnetic potential energy densities almost 
certainly exists in molecular clouds (Zweibel & 
Mckee 1995, and references therein). That is 
3/2po"7VT = By/ 2 /Sit where <7_/vt is the non- 
thermal component of the velocity dispersion, 
which is assumed to arise exclusively from Alfven 
waves. These assumptions lead to the equivalent 
of equation (2) for Bw, that is, 

B w » 0AAv NT nl /2 fiG, (4) 

where Av nt is the non-thermal line width (FWHM) 
and n p is the proton density. Since Zeeman ef- 
fect measurements are generally sensitive to large 
scale fields (on the order of the beam size), they 
are expected to be better indicators of Bs than 
of Bw- Therefore, the formal usefulness of com- 
paring equation (4) for Bw to Zeeman observa- 
tions (Bi os ) is not clear. However, limited obser- 
vational comparisons of Zeeman measurements 
with the Bw predicted from equation (4) show 
that there may also be approximate equiparti- 
tion between the wave's kinetic and static mag- 
netic energy densities in molecular clouds; ie. 
o"tvt ~ l/3fyi ; where v A =Bs 2 /4irp, or equiv- 
alently Bs~Bw ( Myers & Goodman 1988a, 
1988b; Bertoldi k Mckee 1992; Mouschovias & 
Psaltis 1995). In any case it is worth noting that 
Bw must be < Bs to prevent the formation of 
super- Alfvenic shock waves, which are highly dis- 
sipative (McKee & Zweibel 1995). 

We can assess the role of magnetic fields in the 
support of M17 SW from our Hi Zeeman mea- 
surements combined with other measurements of 
the clouds basic properties. Goldsmith, Bergin, 
Sz Lis (1997) estimate that the mass of the M17 
SW is M core = 1.2 x 10 4 M Q (including helium) 
on the basis of their C 18 maps. Using the geo- 
metric mean of the extent of their M17 SW C 18 
maps (1.9 pc x 2.2 pc) ~ 2 pc as the diameter 
of the core region, and their integrated C 18 
column density, we estimate <A r p >« 4 x 10 23 
cm~ 2 and <n p >« 6 x 10 4 cm -3 (assumes fill- 
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ing factor of 1). Although Avnt is often esti- 
mated from the Avfwhm of optically thin lines, 
we believe a more accurate estimate for the in- 
ternal motions of the M17 SW core as a whole 
comes from the velocity dispersion of the cen- 
ter velocities of the clumps themselves. Stutzki 
& Giisten (1990) identified 180, ~ 0.1 pc C 18 
clumps in M17 SW and provide center velocities 
for each. The dispersion of these center veloc- 
ities is 2.7 km s _1 , equivalent to Avnt = 6.3 
km s _1 . From these parameters and equations 
(3) and (4), B s , cr it~ 2000 /xG and B w ^ 620 
/j,G. These estimates indicate that (1) unless the 
static field is four or more times higher than our 
maximum values for Bi os (500 /uG), then M17 
SW as a whole is magnetically supercritical and 
(2) the estimate for By/ is in good agreement 
with our highest values for Bi os . Therefore, the 
evolution of the M17 SW core is not entirely con- 
trolled by Bs and ambipolar diffusion, and in- 
ternal motions must play a role in its support 
and dynamics. Note that these estimates are not 
necessarily valid for individual clumps within the 
core (Bertoldife McKee 1992). 

Another perspective on cloud support can be 
gained from a simplified form of the virial theo- 
rem that simultaneously takes into account mag- 
netic effects from Bs and By/ (but not external 
pressure or rotation). This form is 



W \=M S + ST, 



(5) 



where W is the gravitational energy term, M.$ is 
the magnetic energy associated with Bs, and 3T 
is the energy contribution from all internal mo- 
tions including the magnetic wave energy pro- 
duced by By/ (See, for example, Mckee et al. 
1993 for a more detailed description of these 
terms.). Note that 2T, in the non-magnetic 
virial equation, has been replaced by 3T because 
equipartition between internal motions and Bw 
adds an additional term T to the energetics of 
the cloud. McKee et al. (1993) explain why 
virtually all self-gravitating clouds must be close 
to dynamically critical since their internal pres- 



sures would otherwise be comparable to external 
pressures, contrary to observation. If M17 SW 
is dynamically critical, that is, on the verge of 
collapse, then the virial equation above can be 
rewritten in the form 



B s nBs,crit[l-(3T/\ W|)]V2. 



(6) 



Given the M17 SW cloud parameters estimated 
at the beginning of this section, T/\W\ ~ 0.3, 
so Bs ~ (l/3)Bs,crit ~ 660 fiG, and M$ pa 
1 / 3M core . Given the inevitable errors in estimat- 
ing virial terms, these values are quite uncertain. 
However, they can be regarded as a possible equi- 
librium model for the energetics of the M17 SW 
core. In this model, internal motions in the core 
are sub- virial since T/W < 1/2, and the cloud is 
magnetically supercritical. That is, the cloud is 
supported in part by internal motions, with the 
remainder of the support coming from Bs- More- 
over, the magnitude of Bs in this model is com- 
parable to By/ (620 fiG), so that the magnetic 
support provided by the static and wave compo- 
nents of the field are comparable. This suggests 
that Bs~Byy is a good assumption for M17 SW 
(see above). Although our estimate for Bs us- 
ing this model is quite close to the highest mea- 
sured values of Bi os (pa 500 fiG), this agreement 
must be considered somewhat fortuitous since es- 
timates of Bs are very sensitive to uncertainties 
in cloud parameters and, of course, Bi os is just 
one component of the field. 

Apart from equations (3), (4), and (6), one 
other method has been used to estimate field 
strengths in self-gravitating clouds. Letting 2T= 
a \W\ we obtain the following mass estimate: 



M 



81n2 



AvjfrpR 

aG 



(7) 



where a parametrizes the extent to which a mag- 
netic field is needed to support the cloud, ie. 
a < 1 for clouds in which the internal motions 
are sub-virial. Substituting this estimate for the 
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cloud's mass into equation (4) we obtain 



5. SUMMARY AND CONCLUSIONS 



B\y = 15a 



'1/2 



R 



(8) 



where R is in pc and Avnt is in km s _1 . Note 
that this formulation is similar to those found in 
Myers & Goodman (1988a, 1988b), Mouschovias 
& Psaltis (1995) (do not include a), and Bertoldife 
McKee (1992) (using a). Comparing the mass 
computed from equation (7), M = a _1 8.5 x 
1O 3 M , to the C 18 mass of M17 SW (1.2 x 
10 4 M ), we find a = 0.7 and B w & 740 fiG, 
again, remarkably similar to the peak measured 
value of Bi os . The fact that the field strength pre- 
dicted by equation (8) is relatively well matched 
by B[ os may be another indication that By/~Bs 
for most clouds (see above). 

Although formulations similar to equation (8) 
have often been used to estimate magnetic fields 
in the literature (see above), we consider equa- 
tion (6) to yield a somewhat better field esti- 
mate, at least in so far as comparison with Zee- 
man observations are concerned. This is be- 
cause 1.) if we assume a = 1 ( Myers & Good- 
man 1988a, 1988b), equation (8) is essentially 
the same as equation (4) , both of which estimate 
Bw- There is evidence (this work, and refer- 
ences cited above) that Bw~B$, but this asser- 
tion will require more observational evidence to 
confirm. 2.) In addition, the assumption that 
a=l, excludes entirely the need for support from 
Bs, contrary to our observational Zeeman evi- 
dence, ie. Bi os ~ 500 fiG. 3.) The inclusion of a 
in the equation solves this delima, but creates a 
new one, since the estimate of a only describes 
the degree of under- virialization in the cloud, and 
does not distinguish between support from Bs 
and By/-, while equation (8) is formally an esti- 
mate of Bw alone. For these reasons equation (6) 
may be the best method of estimating the mag- 
netic field based on cloud parameters for the pur- 
pose of comparison with Zeeman measurements 
even though all the estimates found here (ie. eq. 
[4], eq. [6], and eq. [8]) are similar. 



The Hi absorption lines toward M17 show 
complicated profiles with 5 to 8 distinct veloc- 
ity components. In general, the highest Hi 
column densities are concentrated toward the 
Hn region— Ml 7 SW boundary. This effect is 
particularly noticeable in the velocity range of 
the 20 km s _1 Hi component. A lower limit 
to the maximum Hi column density summed 
across the whole Hi velocity range is N^ J -/T S = 
l.OxlO 20 cm~ 2 . Atomic and molecular emission 
temperature measurements suggest that T s lies 
between 50 and 200 K. 

A magnetic field has been detected at the same 
velocity as M17 SW (vlsr ~ 20 km s -1 ) which 
peaks steeply toward the M17 Hn region— M17 
SW interface to values of Bi os = —450 fiG. Analy- 
sis of this H I component's line width, along with 
its velocity and spatial coincidence with molec- 
ular density tracers along the M17 SW molecu- 
lar ridge suggest that it originates in unshocked 
PDR gas. Comparison of the proton density im- 
plied by the observed Bi os and the 20 km s^ 1 H I 
density show that the Hi gas is only a ~ 1% con- 
stituent of the gas along the line of sight toward 
Ml 7. At the same time its distribution along the 
Stutzki et al. (1988) strip scan across the Hn 
region— M17 SW boundary agrees well with that 
of molecular gas. Therefore, we suggest that the 
20 kms H I gas lies in thin, dense shells around the 
numerous molecular clumps which have been ob- 
served in M17 SW, or it must be well mixed with 
fairly dense (1 x 10 4 cm -3 ) interclump molecu- 
lar gas. In addition, the region of maximum Bi os 
for this component agrees spatially with a min- 
imum in 100 fim polarimetry which samples the 
transverse magnetic field. 

Another significant Bi os was observed in the 
blended component at 11—17 km s _1 with values 
reaching ~ + 550 /xG. The morphology of high 
Nhi for this component is quite different from 
that at 20 km s _1 , with the greatest concentra- 
tions found further east and toward the north- 
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ern parts of the source. From these components' 
morphology, line width, and evidence that there 
is a significant amount of obscuring material on 
the front side of the M17 Hn region, we suggest 
that the 11 — 17 km s _1 components originate in 
shocked gas that is streaming toward us along 
the line of sight. Estimates of the shock speeds 
which can be produced by the M17 Hn region 
agree well with these components' blueshift (2—8 
km s _1 ) with respect to the rest velocity of the 
region. Like the 20 km s _1 Hi component, the 
11—17 km s _1 Hi gas can only comprise ~ 1% 
of the line of sight gas, and may be confined to 
photodissociated Hi shells surrounding molecu- 
lar clumps close to the Hn regions front side. 

Using various virial arguments we have esti- 
mated that M17 SW is sub- virial (T /W< 1/2) 
and magnetically supercritical (M core > M$). In 
our model ~ 1/3 of M17 SW's total support is 
provided by its static magnetic field B$ (~ 660 
fj-G), with the rest of the support arising from in- 
ternal motions including support from the wave 
component of the field B\y. Estimates of the 
static (Bs) and wave (Byy) components of the 
magnetic field indicate that they are of approx- 
imately the same magnitude (~ 600 — 700 fiG) 
and agree well with our highest values of Bi os 
(» 500^G). 

Subsequent analysis of higher resolution VLA 
Hi Zeeman effect data will give better sensitiv- 
ity to Bi os if the field is tangled, and allow us 
to better differentiate separate velocity compo- 
nents. In addition, observations of the Zeeman 
effect in OH may also help untangle the magnetic 
field data at this interface since it tends to trace 
somewhat higher density gas. 
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Table 1 

Observational Parameters of 21 cm HI Zeeman Observation of M17 



Value 



1420 MHz 
1996 June 2 
8 hr 
30' 

61.0" x 44.5" 
18 h 17 m 35.0 s , -16°14'00" 
256 

781.25 kHz (165.05 km s -1 ) 
-67.5 to +97.5 km s -1 
3.052 kHz (0.64 km s -1 ) 
40 mJy beam -1 
80 mJy beam -1 
1 mJy beam -1 = 0.25 K 
T= 0.6 pc 
18 h 17 m 32.5 s , -16°13'06" 



Parameter 



Frequency 

Observing date 

Total observing time 

Primary beam HPBW 

Synthesized beam HPBW 

Phase and pointing center(B1950) 

Frequency channels per polarization 

Total bandwidth 

Velocity coverage 

Frequency resolution 

Rms noise in line channels 

Rms noise in continuum 

S v to T b 

Angular to linear scale"" 

Position of 21 cm continuum peak 



a Assumes distance of 2.2 kpc. 
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